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ABSTRACT 



Finite element, elastic and plastic analyses of various 
aluminum panels, containing holes and notches, were conducted 
for comparison with photoelastic experimental results. A 
FORTRAN IV program, ADINA (Automatic Dynamic Incremental 
Nonlinear Analysis) , was used for both linear and nonlinear 
analyses. Mesh refinements were used for each panel and the 
monotonically convergent results were extrapolated using 
Richardson's method. Stresses were locally smoothed from 
the Gauss integration points to the nodal points. Eight 
noded, isoparametric elements were used throughout. Modi- 
fication to an ADINA preprocessor program, also coded in 
FORTRAN IV, was made for use with a VERSATEC plotter. 

Comparisons were made to the elastic, analytic series 
solution by Howland for a circular hole in a finite strip. 

The finite element results varied by less than one percent 
from Howland's solution. Handbook values for the elastic 
stress concentration factors of the geometries investigated 
differ from finite element results by less than one percent 
in all cases. The photoelastic works of Frocht were also 
compared where applicable. Stresses in the plastic range 
obtained from slip-line theory for a rigid-perfectly-plastic 
material show excellent correlation to a finite element 
analysis of such a material. Comparisons to elastic and 
plastic experimental data were made for the panels 
analyzed and show good correlation to finite element results. 
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I. 



INTRODUCTION 



Development of on-board fatigue monitoring systems 
for Naval aircraft have made it possible to record exten- 
sive structural loading data in flight. The strain gages 
used in such a system must be located away from stress 
concentration areas to prevent their fatigue; however, 
these areas are of the greatest interest in analyzing 
and predicting fatigue life of the structure. Understand- 
ing the relationship between nominal, far-field stresses 
and local stresses in critical areas thus becomes vitally 
important. Recent experimental investigations into the 
effect of uniform, far-field loads on stress concentration 
areas have been made at the Naval Postgraduate School 
(NPS) using photoelastic techniques [Refs. 1, 2 and 3]. 
These experiments involved loading 7075-T6 aluminum into 
both the elastic and plastic regions, as well as measure- 
ments of residual stresses resulting from plastic yielding. 

Finite element analyses (FEA) of the aluminum panels 
used in the experiments of Stenstrom [Ref. 1] were con- 
ducted. The panels used in the experiments of Engle 
[Ref. 2] and Stuart [Ref. 3] have similar geometry. The 
finite element programs available at NPS were surveyed 
and ADINA [Ref. 4] was chosen for its proven ability to 
produce the nonlinear analyses required for plastic 
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yielding of aluminum. To provide increased accuracy, 
each panel was modeled using two meshes. The results 
obtained for the coarse and fine meshes were extrapolated 
to a final result using the Richardson extrapolation 
technique [Refs. 5 and 6]. 

Along with ADINA, a preprocessor program, PSAPl [Ref. 7], 
was used to verify mesh connectivity prior to analysis by 
ADINA. PSAPl provides a graphical output of the finite 
element mesh and was coded for the CALCOMP plotter installed 
at NPS prior to 1978. For this thesis PSAPl was adapted 
for use with the VERSATEC plotter now installed at NPS. 

The stress-strain material properties of the 7075-T6 
aluminum actually used to make the panels had to be estab- 
lished to provide an accurate material model for use with 
ADINA. Material testing was conducted to establish the 
Young's Modulus (E) , Poisson's Ratio (v) , yield stress (o^) , 
strain hardening modulus (E^) and the Ramberg-Osgood co- 
efficients 6 and n. 

Comparisons were made to other works, in addition to 
the experiments conducted at NPS. The initial analysis 
involved a comparison of FEA to the results of Howland 
[Ref. 8], for a circular hole in a finite strip, to vali- 
date the methods used. A comparison of FEA to plane 
stress, slip-line theory, for rigid-perfectly-plastic 
material was also included as a validation for the plastic 
analyses. 
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II. MATERIAL PROPERTY TESTING OF 7075-T6 ALUMINUM 



The elastic and plastic material properties of the 
aliominum panels were established by tensile tests of uni- 
axial specimens made from the same mill run. The speci- 
mens were manufactured and tested according to current 
ASTM standards [Ref. 9]. MICRO-MEASUREMENTS, EA-13-125AD- 
120, precision strain gages with a temperature compensated 
bridge circuit were used on all specimens. Transverse 
gage sensitivity errors were corrected according to the 
manufacturer's recommendations [Ref. 10]. Critical cross- 
section measurements were made with a micrometer. 

A. TESTS FOR AXIAL LOADING 

Initial tests of the two-gaged specimen. Fig. 1, in 
the MTS testing machine indicated a significant bending 
moment was being produced by the 30,000 lb GRIFF grips. 

To investigate this problem further, tests were conducted 
on both the MTS and RIEHLE test machines with a five- 
gaged specimen shown in Fig. 2. The results of these 
axial loading tests, shown in Table I, verified that the 
GRIFF grips on the MTS test machine do not give axial 
loading. An inspection of the gripped region on the speci- 
men indicated that the jaws of the grip were not applying 
a uniformly distributed force and thereby induced a bending 
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moment by off-axis loading as shown in Fig. 2. The grips 
on the RIEHLE test machine gripped evenly and a uniform 
strain distribution resulted as seen in Table I. 

B. CHARACTERISTICS OF 7075-T6 ALUMINUM PANELS 

The following characteristic properties were determined 
from the four specimens tested. 

1 . Young's Modulus (E) 

Tests were conducted using the specimen shown in 
Fig. 3 on the MTS test machine with 10,000 lb INSTRON grips, 
which gripped the specimen evenly. The results of testing 
three specimens are shown in Tables II to IV. Linear 
regression in the elastic range of all the test results 

g 

determined a Young's Modulus of 10.12 x 10 psi, with a 
correlation coefficient of 0.9996. 

2 . Poisson's Ratio (v) 

Tests were conducted using the specimen shown in 
Fig. 1 on the RIEHLE test machine with 10,000 lb RIEHLE 
grips. The results are tabulated in Table V. Linear 
regression of these results in the elastic region deter- 
mined Poisson's Ratio to be 0.3256 with a correlation 
coefficient of 0.99996. 

3 . Yield Stress and Strain Hardening Modulus 
These values, required for ADINA's bi-linear 

material model, were detemined graphically using the 
data from Tables III and IV. Plastic region data in 
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Table II is not reliable because of excessive creep en- 
countered during that test. 

0.2% offset yield stress, = 76,000 psi 
strain hardening modulus, = 566,000 psi 
The graphical fit of these values to the test data can be 
seen in Fig. 4. 

4 . Ramberg-Osgood Coefficients 

The Ramberg-Osgood equation for elastic-plastic 
stress-strain characterization is given by: 



= 1 ^ 



B(|)" 



( 1 ) 



where : 

e = strain 
a = stress 
E = Young ' s modulus . 

The B and n coefficients were determined graphically from 
the data of Table IV, by the method given by Rivello [Ref. 11] . 
The data in Table IV gave the following values which are 
the best fit to the combined test data 
B = 1. 479 X lo"^^ 
n = 21.58 

The graphical fit of these values, in Eq. (1) , with the 
test data is also shown in Fig. 4. 
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III. MODIFICATION TO GRAPHICAL PREPROCESSOR 



The use of a graphical preprocessor program, such as 
PSAPl, is vital in detecting mesh errors that may other- 
wise go unnoticed. Establishing the proper node locations 
and element geometry prior to analysis for a complex code 
such as ADINA is of utmost importance. 

A. PSAPl MODIFICATIONS 

The program PSAPl was originally coded in FORTRAN IV 
for use on the NPS IBM 360/370 installation with the 
CALCOMP Model 765 drum plotter. The CALCOIiP system once 
installed at NPS used the +Y axis as the unlimited plotting 
axis, see Fig. 5. The entire plotting logic in PSAPl uses 
this orientation of axes to allow multiple plots in a con- 
tinuous strip. With the VERSATEC Model 8222A electrostatic 
plotter now installed at NPS, the +X axis becomes the un- 
limited plotting axis, shown in Fig. 5. To avoid an ex- 
tensive recoding of PSAPl for use with the VERSATEC plotter 
a simple coordinate transformation of the plot was made in 
a limited number of short subroutines. First, all installa 
tion dependent plotting calls used in PSAPl were identified 
These involved seven plotter functions for which new sub- 
routines were coded - 
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Function 



New Subroutine 



Initialize Plotter 



CALCMP 



Move Plotter Pen 



CALPLT 



Letter on Plot 



NOTATE 



Number on Plot 



CALNUM 



Determine Current Pen Location 



CALWH 



Draw a Line 



CALINE 



Stop Plotter 



PSTOP 



The subroutines listed above merely rotates the plot 
to coincide with the VERSATEC axis orientation and retain 
all features originally in PSAPl. Since all plotter 
hardware code is now isolated in these seven subroutines, 
future adaptations to other plotting systems is simplified. 
To provide documentation of this update to PSAPl, a complete 
listing of the new program is provided in Appendix D. 

B. USE OF PSAPl 

Previous use of PSAPl on the IBM 360 system necessitated 
use of a load module since PSAPl took over one minute of 
CPU time to compile. With the new IBM 3033 system 
compilation requires eight seconds; however, use of a load 
module or disk stored source code is still recommended 
since PSAPl contains roughly 2,500 lines of code. Appendix 
A contains sample JCL to use PSAPl on the IBM 3033 MVS 
system. With minimal effort PSAPl could also be set up 
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for use on the IBM VM/370 system. The user's manual for 
PSAPl is in Ref. 7. In addition to a mesh plot, PSAPl 
provides a listing of the node coordinates, element 



connectivity and several key input values 
tion of ADINA. This information provides 
of the input data. 



used in execu- 
a useful check 
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IV. FINITE ELEMENT ANALYSIS (FEA) 



A. DESCRIPTION OF MESHES USED 

1 . Length to Width Ratios and Boundary Conditions 
Initial finite element models of specimens had 
length to width ratios near one, as used by Garske [Ref. 12] 
for his FEA, but they did not provide the desired uniform 
distribution at the loading boundary. Specimen length to 
width ratios of 3-5 were used by Armen, Pifko and Levine 
[Ref. 13] in their FEA and by Stenstrom [Ref. 1] in his 
photoelastic experiments. The criteria established to 
determine uniform boundary stress distribution was uni- 
formity in nodal displacements along the loaded edge as 
discussed by Segerlind [Ref. 14]. In the models used for 
FEA in this thesis, nodal displacements were uniform to 
within 0.1%, and the resulting stress distribution was 
uniform axially to within 0.1% at the panel ends. In all 
cases two-dimensional, eight noded, isoparametric elements 
were used. These higher order elements cannot be loaded 
in an "intuitive" manner as discussed by Zienkiewicz 
[Ref. 15, p. 223]. Figure 6 shows the nodal loading 
required to obtain a uniform surface load. 
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2. Element Meshes 



Two meshes were developed for each panel analyzed. 

A reasonable effort was made to keep element corner angles 
as close to 90° as possible to reduce the effect of element 
distortion discussed by Hopkins and Gifford [Ref. 16] . 

All meshes modeled a quarter of the actual panel by using 
the two axes of symmetry as is common practice in FEA. 

The step from course to fine element meshes was made so 
that each element in the course mesh was subdivided into 
four smaller elements of the same type. Such a mesh sub- 
division can be expected to give monotonic convergence of 
results, Cook [Ref. 17], and allow extrapolation to results 
of an infinitely fine mesh. Figures 7 through 13 illus- 
trate the element meshes used in this analysis as plotted 
by PS API. 

B. COMPUTATIONAL PROCEDURES 
1. Using ADINA 

Once the mesh has been developed, input data is 
prepared in accordance with the ADINA user's manual [Ref. 4]. 
This same set of data is then used as input for PSAPl to 
check for errors and provide a graphical display of the 
element mesh. After preprocessing by PSAPl, the data is 
entered into ADINA for analysis. In the case of linear 
analysis, two types of stress output may be specified, 
nodal point or Gauss integration point. Nodal point output 
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can be computed for up to eight node point stresses for 
each element. Since 2x2 Gauss integration was used, four 
Gauss point stresses were computed for each element. The 
2x2 Gauss integration is recognized as the most efficient 
integration order for this type of analysis [Ref. 15, p. 284]. 
The linear analysis used an isotropic linear elastic material 
model (MODEL "1" in Ref. 4) which required input of E and 
V material properties. The nonlinear analysis allows only 
Gauss point stress outputs and uses a bilinear elastic- 
plastic material model, with von Mises yield condition and 
isotropic strain hardening (MODEL ”8" in Ref. 4) . 

For static analyses ADINA uses a time function 
method to apply loads in steps. Linear analysis loading 
was accomplished in a single step to a nominal value of 
3,000 lbs load. Nonlinear analysis loads were applied in 
ten steps to a maximum value , matching the experimental 
loads, and then unloaded to zero in ten steps to obtain 
residual stresses. The stress output from ADINA is a 
listing of nodal or Gauss point stresses for each element. 
Since the only area of interest in this analysis was the 
distribution of stresses along the reduced cross-section, 
no large post-processing program was developed or used. 

All final computations using ADINA output data were accom- 
plished on a HEWLETT-PACKARD 9830A calculator, using short 
programs coded in BASIC. If more extensive stress 
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distribution information were desired, some form of auto- 
mated post-processing would be necessary to reduce the 
computational workload. At a minimum, nodal stress outputs 
by ADINA must be averaged to obtain unique values of stress 
at nodes shared by more than one element. 

2 . Richardson Extrapolation 

The use of course and fine meshes allows extrapola- 
tion to an infinitely fine mesh as discussed earlier. 
Richardson extrapolation [Ref. 5] was used in this analysis 
where : 

a = — b — £ £ c — (2) 

extrap h ^ - h ^ 

“f ^c 



where 



a 

a 

a 



extrap 

c 

F 




= extrapolated solution 
= solution obtained with h=h 

c 

= solution obtained with h=h„ 

r 

= linear dimension of course element 
= linear dimension of fine element 



m =2 (for this analysis) 

The exponent m is determined by the order of the discre- 
tization error O(h^). Since h represents the length of an 

2 

element the element area is represented h . In a two 

dimensional problem such as this O(h^) is of the order of 
2 

h , the area of an element. In the mesh refinement scheme 
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Equation (2) can be rewritten 



used h_ 

r 



1 , 

2 ‘"c 




thus 



a 



extrap 



a 



extrap 



2 % 2 

“c<iT' - “f‘5T) 

c c 

hp 2 h 2 

<ir> - <iT> 

c c 




(3) 



(4) 



Equation (4) then becomes the relation to obtain extrapolated 
stresses from coarse and fine mesh results in a two dimen- 
sional analysis. Better extrapolations can be obtained by 
using three or more refined meshes, but the computational 
effort increases significantly. 

3 . Optimal Stress Locations and Local Smoothing 

It is generally accepted that the most accurate 
sampling points for stresses are the Gauss integration 
points within the element [Ref. 15, p. 281, and Ref 18]. 

In this analysis, the nodal points are of the greatest 
interest; thus a technique of local smoothing must be 
applied to the integration point stresses to obtain nodal 
stresses as reported by Hinton and Campbell [Ref. 19] . 

The formulation of this local smoothing technique for ADINA 
elements is developed in Appendix B. The nodal values 
obtained must then be averaged if shared by two or more 
elements. 
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4 . Computational Times 



Because of the extremely large size of ADINA (about 
17,000 lines of code in the NFS version) and the out of core 
solver, it does not adapt well to time sharing systems. 

Using the IBM 360 system at NFS, ADINA required 31 user 
defined overlays to create a manageable load module in 
about 30 minutes of CFU time. With the new IBM 3033 MVS 
system at NFS, ADINA is compiled without overlays in about 
one minute. When using a load module, the program execu- 
tion took less than 2 minutes CFU time on the IBM 3033. 

In addition to ADINA, the preprocessor (FSAFl) 
and post-processing techniques involve considerable time 
and effort. Figure 14 is a flow chart of the computational 
procedure used in this analysis. An example of the JCL 
to use ADINA in load module form on the IBM 3033 MVS 
system v;ith use of the mass storage facilities is shown 
in Appendix C. 
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V. RESULTS OF ANALYSIS 



A. CIRCULAR HOLES IN LINEAR MATERIAL 

The FEA results for a circular hole in a finite width 
strip were used to validate the elastic computational pro- 
cedure discussed earlier. The results of Howland [Ref. 8] 
were compared to both the Gauss point smoothed results and 
the nodal output results in Fig. 15. The stress concentra- 
tion factor a/a is referenced to the far-field stress. 

The smoothed results give the best match to the results 
of Howland at the edge of the hole, and the only signifi- 
cant variation between the two FEA methods occurs within 
the first 0.25 inches from the edge. In order to obtain 
the 0.25 inch stress value for the coarse mesh, in the 
Gauss point smoothed result, a midside node value had to 
be obtained by the averaging method discussed in Appendix 

B. The linear distribution of smoothed stresses along the 
sides of the element, [Ref. 19] , appears to produce a 
less accurate result in this area of extreme stress 
gradient, when compared to ADINA's nodal output result. 

This tendency was noted in all cases; however, the peak 
stress values from smoothed results consistently gave 
better correlation with other investigators [Ref. 20]. 

A circular hole with 1=0.25 was also analyzed and 
compared to the experimental results of Stenstrom [Ref. 1] 
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along with an interpolation of Howland's results. The Og 
experimental data correlates well with the FEA results; 
however, the experimental data shows significant varia- 
tion between 0.125 and 0.375 inches from the edge of the 
hole, as seen in Fig. 16. 

B. OPPOSITE U NOTCHES IN LINEAR MATERIAL 

The results for linear analysis are presented in non- 

dimensional stress concentration form; however, the 

normalizing stress changes. For U notches is the 

stress concentration factor referenced to the theoretical, 

nominal stress (a ) in the reduced cross-section where 

n 

= Load/ Area of Reduced Cross-Section. 

1. Shallow Notch Panel 

The FEA results plotted with the experimental 
data of Stenstrom are shown in Fig. 17. Once again both 
FEA results are shown and the variation for the two methods 
occur within 0.25 inches from the notch edge; however, 
there was less variation than was seen in the circular 
hole analyses. 

For this panel the experimental data appear to be 
uniformly below the FEA results for Cg. The data shows 
significant variation at the 0.625 inch point but follows 
the proper trend within 0.5 inches from the notch edge. 
According to data collected by Peterson [Ref. 20] , this 
notch geometry should yield a maximum = 2.74. The 
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Gauss point smoothed results matched this value exactly. 
Stuart [Ref. 3] reported a = 2.69 for this same notch 
geometry, with a standard deviation of 0.187 in the 14 
samples he measured by photoelastic methods. Early photo- 
elastic work by Frocht [Ref. 21] detemined = 2.7 for 
this notch geometry. The FEA results appear to be in good 
agreement with other investigators^ for this notch geometry. 

2 . Deep Notch Panel 

Results of the FEA for a deep notch were plotted 
with Stenstrom's experimental data in Fig. 18. The results 
of the two FEA methods again diverge within 0.5 inches from 
the notch edge. FEA stress values at the 0.25 inch point 
have spread farther apart in this case since the stress 
gradient is very severe at that point. The experimental 
data correlates well for both and a ; however, the maxi- 
mum experimental Og is considerably lower with a = 3.83. 
Results reported by Stuart for this notch geometry was a 
= 4.05 with a standard deviation of 0 . 219 for 14 speci- 
mens measured by photoelastic methods. Frocht [Ref. 22] 
reported a photoelastic = 3.9 for this notch geometry, 
but concluded that the result was 5-10% low, giving a 
corrected range of from 4.1 to 4.3. The Gauss smoothed 
FEA result gave a K,p = 4.24 which compares well to an 
empirical relation given by Peterson [Ref. 20] for 
r/d < 0.25. 
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(0.78 + 2.243/t7r) 



0.993 + 0.18 ^ 



2t 2 2t 3 

1.06(^) + 1.71(^) 



where 



t = notch depth (3.9375) 
r = notch radius (0.625) 
d = minimum width (15.625) 

D = maximum width (23.5) 

Inserting the values above into Eq. (5) gives a = 4.26, 
which is 1/2% above the FEA result. Other FEA results 
reported by Armen, Pifko and Levin [Ref. 13] and Griffis 
[Ref. 23], using linear strain triangle (LST) elements, 
produced K,^ values within 5% of those produced by use of 
Eq. (5) . The rectangular elements used in this analysis 
are known to give better results than LST elements as 
noted by Clough [Ref. 24] . It is clear that the FEA 
results obtained for this notch are in good agreement 
with other works. 



C. OPPOSITE U NOTCHES IN NONLINEAR MATERIAL 

The analysis for loading into the plastic region of the 
7075-T6 aluminum was made using the bilinear material 
model discussed earlier. The loads used were selected to 
match those used in the experiments of Stenstrom; thus 
allowing direct comparison. The strains obtained in those 
experiments were used to solve for stresses by use of the 
Prandtl-Reuss plastic flow equations. 
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1. Shallow Notch Panel 



The results of FEA for the three load cases, 60,000, 
65,000 and 70,000 lbs are presented along with the experi- 
mental results in Figs. 19 through 21. The Cg results 
compare well although no trend for peak Og stress away from 
the notch edge is shown in the experimental data. In all 
cases the FEA determined the peak Og stress to occur near 
the yield boundary, and the gradient of the Cg stress to 
fall off dramatically in the plastic zone. This character- 
istic behavior of the Cg stress was reported by other 
investigators [Refs. 13 and 23] using FEA on 2024-T3 alxominum. 
Plane elastic-plastic stress distributions reported by 
Frocht [Ref. 25] show similar trends. The experimental 
data also shows a marked change in the gradient of Cg stress 
within the plastic region. The growth of this plastic 
region is approximated using the FEA results for this notch 
in Figs. 22 through 24. 

Experimental data for the stress distribution 
matches the FEA results closely except at the notch edge 
where the measured does not go to zero as it should. 

The characteristic peak value of near the plastic 
boundary as seen in the FEA results is also shown by the 
test data. 

The FEA residual stress computed upon unloading 
from the three load cases are shown in Figs. 25 through 30. 
The characteristic distributions of the residual stress 
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agrees with those reported by others [Refs. 25, 26 and 27]. 

The experimental residual stress distributions reported by 
Stenstrom show similar trends but significant variations 
when compared to the FEA results . 

2 . Deep Notch Panel 

Three load cases were computed to plastic loading 

levels; however, only limited experimental results are 

available for this notch as seen in Fig. 31. The 30,000 lb 

load is just at the onset of yield in the notch root area. 

Limited residual experimental data [Ref. 3] was available 

for comparison in Figs. 32 and 33 which are plots of the 

residual and a stress distributions as a result of the 
9 r 

three loading cases. The 30,000 lb load has caused yielding 
in a small region at the root of the notch as seen in Fig. 34. 
Figures 35 and 36 illustrate that the plastic zone does not 
grow to the extent it did in the shallow notch. The stress 
gradients are very severe in the deep notch and as a result 
the 0.25 inch sampling points used in the FEA may be useful 
in only showing gross trends close to the notch edge. The 
trends appear to be much the same as in the shallow notch, 
with the peak o stresses occurring near the yield 

boundary. The yield boundaries shown in Fig. 31 are 
approximations based on qualitative analysis of the finite 

element results. The a, experimental data for the 30,000 

0 

lb load case correlates well with the peak stress, again 
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appearing low as it did in the linear analysis. The residual 

0Q stress distribution in Fig. 32 follow much the same trends 

as seen in the shallow notch but with much higher gradients 

within the first 0.5 inches from the notch. The FEA results 

for the residual stresses shown in Fig. 33 indicate a 

limitation in the element size used since the a value of 

r 

the notch edge does not return to zero as it should. Because 
of this problem the data may be questionable for showing 
proper trends in the first 0.5 inch from the notch edge. 

3 . Rigid-Perfectly Plastic Panel 

A stress distribution for the theoretical material 
used in slip-line theory was desired. By using ADINA's 
bilinear material model such a material could be approxi- 
mated reasonably well. For this analysis a Young's Modulus 
2 6 

(E) of 10 psi was used to model perfect rigidity. The 
strain hardening modulus (E^) was set to zero to model 
perfect plasticity. Poisson's ratio (v) was 0.4999999, 

2 6 

as close to 0.5 as the computer would allow. The E = 10 
also represents a computer limit in approaching an infinitely 
large E. Figure 37 illustrates the results obtained and 
compares them to a slip-line solution. The results are 
normalized to the arbitrary 73,000 psi yield stress used 

in the analysis. The g values obtained agree exactly 

0 

with slip-line theory. Conversely the results do not 
reflect the same values as slip-line theory, but do show 
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a similar trend. The growth of the plastic zone obtained 
is shown in Figs. 38 through 40. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



The results obtained from the FEA have proven useful 
in determining the validity of experimental data gathered 
by photoelastic techniques. The FEA results varied by 
less than 1% when compared to published analytical results 
and handbook values. Experimental data from Stenstrom's 
photoelastic work correlated well with the FEA results. 

The primary exception was the residual stress experimental 
values, which varied significantly from the FEA results. 

The variation may be in transforming the residual strains 
measured photoelastically into stresses for comparison 
with the FEA results, since ADINA only provides a stress 
output. Limitations of ADINA' s bilinear material model, 
initially considered severe, do not appear to have hampered 
this investigation. A possible exception is the residual 
analyses where the transition region from elastic to 
plastic strains becomes especially important. The Gauss 
point smoothed results gave the best correlations at the 
edge singularities in all cases; however, due to the 
limitations noted at 0.25 inches from the edge, the 
results at that specific point may not be as accurate for 
this method. The nodal output results gave consistently 
higher stress values at the edge singularities. Because 



35 



of the severe stress gradients near the deep notch analyzed 
the use of a finer element mesh near the notch would probably 
produce better results. 

The effort involved in developing two meshes such as 
those used in this thesis is considerable. An automatic 
mesh generation capability would reduce the workload and 
allow experimentation with several types of element meshes. 

ADINA proved to be a useful and powerful program, as 
expected, but something simpler and less awkward to use may 
be all that is required for two dimensional analysis. Such 
a system is already in use at NPS but does not offer non- 
linear capabilities. If use of ADINA is to be continued in 
this type of investigation, a post-processing program should 
be adapted. There are programs available to post-process 
ADINA data at NPS [Refs. 28 and 29] but they would require 
modifications to work with two dimensional analyses and 
the VERSATEC plotter. 

Standardized material property testing would ease the 
inevitable task of obtaining basic material properties for 
use in analysis or experiments. Some form of automatic 
data collection with use of the MTS testing machine would 
allow testing of a larger sample population and provide 
statistically more accurate information. 
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FIGURE 1 



2 GAGE SPECIMEN 
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FIGURE 2 
5 GAGE SPECIMEN 
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FIGURE 3 
1 GAGE SPECIMEN 
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FIGUPE 5 

CALCOKP AND VEPSATEC PLOTTER AXES 
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FIGURE 6 

NODAL LOADING DIAGRAM 
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FIGURE 7 

COURSE MESH FOR CIRCULAR HOLES 
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FIGURE 8 

FINE MESH FOR CIRCULAR HOLES 
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FIGURE 9. COURSE MESE FOR SHALLOW NOTCH 
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FIGURE 13 

EXATIPLE OF COMPLETE PANEL MESHES 
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FIGURE 22 

SHALLOW NOTCH 60,000 LB LOAD PLASTIC ZONE 
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SHALLOW NOTCH 65,000 LB LOAD PLASTIC ZONE 
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FIGURE 24 

SHALLOW NOTCH 70,000 LB LOAD PLASTIC ZONE 
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FIGURE 34 

DEEP NOTCH 30,000 LB LOAD PLASTIC ZONE 
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FIGURE 35 

DEEP NOTCH 35,000 LB LOAD PLASTIC ZONE 
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FIGURE 36 

DEEP NOTCH 40,000 LB LOAD PLASTIC ZONE 
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FIGURE 38 

RIGID-PERFECTLY-PLASTIC INITIAL PLASTIC ZONE 
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FIGURE 39 

RIGID-PERFECTLY-PLASTIC INTERMEDIATE PLASTIC ZONE 
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FIGURE 40 

RIGID-PERFECTLY-PLASTIC FINAL PLASTIC ZONE 
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TABLE' I. 



MTS AND REIHLE 5 GAGE TEST RESULTS 



““6 

All Strains are 10 in/in 



MTS Test Machine 



LOAD 

lbs 



STRAINS 



1000 


-3 


421 


814 


1235 


1638 


2000 


338 


1003 


1619 


2279 


2903 


3000 


960 


1714 


2431 


3186 


3908 



REIHLE Test Machine 



LOAD 

lbs 



STRAINS 



1000 


800 


763 


755 


741 


726 


2000 


1564 


1540 


1546 


1543 


15 35 


3000 


2370 


2334 


2363 


2377 


2388 
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TABLE II . MTS SPECIMEN A TEST RESULTS 

2 

Cross-section = 0.03975 in 



Load , lbs . 


“6 

strain , e , 10 in/in 


25b 


615 


503 


1204 


750 


1801 


1005 


2423 


1255 


3042 


1505 


3665 


1778 


4352 


2003 


4925 


• j • J cr 


C C" 

•-I • J C- \ 


2508 


6355 


2755 


7365 


2905 


8 2 3 0 


2 9 3 4 


9045 


3037 


1 0 1 50 
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TABLE III. MTS SPECIMEN B TEST RESULTS 



Cross-section = 0.03975 in 



Load, Ibs- 


~6 

Strain, e , 10 in/in 


500 


H ■ iz: I .i 


ir |-i 


1 8S0 


1 j y 


1 0 0 0 


2450 


1250 


3030 


1500 


3700 


1750 


4320 






2 0 0 0 


'4 9 5 0 


2250 


J b 2' y 


2500 


6330 


2750 


r'340 


2900 


3200 






2000 


9250 


3 0 5 0 


1 0200 


3 1 0 0 


1 1550 


3125 


1 2800 
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TABLE IV. MTS SPECIMEN C TEST RESULTS 

2 

Cross-section = 0.03975 in 



Load, lbs. 


—6 

Strain, £,10 in/in 


272 


b 6 2 


503 


122U 


7 b 2 


1850 


1 0 0 8 


2451 


1231 


3 0 0 5 


150b 




1755 


431 3 


2000 


•4 '2 ‘‘t i J 


■ j “i cr irj 


J b0 i’ 


2503 


b' o 2 


2750 


7 2 ^5 i-3 


2900 


:3 1 2 0 


300 0 


9 2 0 0 


3 0 b 0 


i 

cn 


3 1 0 0 


I 1 500 


3125 


.!. 2 5 0 0 
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TABLE V 



REIHLE SPECi:4EN TEST RESULTS 



Cross-section 



0.12 



. 2 
in 



LOAD 


STRAIN, Si 


STRAIN, £2 


lbs. 


10 ^in/in 


10 ^ in/in 


5 0 0 


345 


1 1 ’5 


1 000 


720 


-2 '^10 


1500 


1 1 1 4 


-384 


2000 


1501 


-492 


2 5 0 0 


1895 


" ' *3 1 '4 


I-: 0 0 0 


3291 


4 3' 


3500 


2705 


..... '7 1 


4000 


3095 


'P 


4500 


3515 


i i'j 


5 0 0 0 


3912 


-1255 


5500 


4332 


- 1 : 4 


60010 


4 750 


1 5 3 i-l 


6500 


cr ■» ■ 1 “ 
J i '1' ' ’j 


-i.-:.4? 


7000 


cr cr c_- 


..-17:34 


7500 


8075 


- L'-523 


8000 


6849 


-2 1 03 


3 5 01 0 


CT 

r 2 0-J 


1-2 1 _• 


800101 


6 6 3 


- 2 3 1: 8 


9 5 01 01 


12245 


-4 4 35 



81 



Ti^BLE VI. 



TABLE VII. 



X = 0.2 HOWLAND DATA 



DISTANCE 
FROM HOLE 
in . 




0 . 0 


. 1 4 


0. 3 


1 . 57 


1 . 0 


1 . 26 


1,3 


1 „ 1 6 


2, 0 


1 . 1 1 


2,3 


1 . 07 


:3 , 0 


1 . 05 


•”i cr 

,"i , .j 


1.01 


4, 0 


0. 97 



X = 0.25 HOWLAND DATA 

DISTANCE ^ 

FROM HOLE 9/ 0 ^ 

in. 



0 . 000 


■J . c 6 


0. 422 


1 . 75 


0. 828 


1 . 30 


1 . 234 


1 . 22 


1 641 


1 . 1 .? 


2. 047 


1 . 07 


2.458 


1 . 04 


2. 859 


0. 97 


3 . 063 


0 . '? 5 
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TABLE VIII 



X = 0.2 FEA RESULTS - NODAL OUTPUT 



DISTANCE 
FROM HOLE, in. 


a = 5000 

oo 

'^ 0 / psi 


psi 




f psi 


0 . 00 


15038.2 






351 , 7 


0.25 


10259,7 






2375, 8 


0.50 


"7 cr ■“:» 

1 r J .. t 






1 1 ' c' -9 s b 


0 . 75 


O •“! -j 01 






1673, 1 


1.00 


6307.4 






1 334 , 7 


1.25 


5995.4 






1087, 2 


1 . 50 


5797, 7 






874. 4 


1 . 75 


5b 33 . 9 






700 . 5 


2 . 0 0 


5509 . 3 






5 5 6 , i3 


TABLE IX. 


X =0.2 FEA 


RESULTS 


- GAUSS OUTPUT 




= 5000 


psi 






DISTANCE 
FROM HOLE, in. 


/Psi 






,psi 


0 . 0 0 


15595,4 






128.9 


0. 25 


9 4 8 c . 6 






4 1 0 4 . 1 


0 . 5 0 


7953.3 






1 326 . 6 


0 . 75 


6 7 3' 3! . 7 






1 650 . 4 


1 . 0 0 


■6316, 1 






1 348, 7 


1 . 25 


5975.6 






1088, 8 


1 . 50 


5795.8 






880,. 1 


1 . 75 


5623 . 8 






695 . 7 


2 . 0 0 


5499.9 






554.9 
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TABLE X 



X = 0.25 FEA BESULTS - NODAL OUTPUT 



DISTANCE 



a = 4500 psi 



FROM HOLE, in. 


Oq fPsi 


Oj. /Psi 


0 . 00 


14745.5 


■h S 9 . 3 


0.25 


9 4 4 7 - 0 


2 1 ;5 1 . 1 


0 . 5 0 


7 ii 3 2 . 6 


1 599 . 3 


0. 75 


6304 . 7 


1473.5 


1 . 00 


5771 . 7 


1 1 39 . 6 


1 . 25 


5471 . 2 


391 . 6 


1 . 50 


5266. 9 


67 1 . 0 


1 . 75 


5091 . 7 


4 y 9 . 1 


2 . 00 


4940. 2 


354.5 



TABLE XI. X = 0.25 FEA RESULTS - GAUSS OUTPUT 

a = 4500 psi 

DISTANCE 
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TABLE XII 



SHALLOW NOTCH FEA LINEAR RESULTS - NODAL 





a = 4225 psi 
n ^ 




DISTANCE 






FPOM NOTCH, in. 


,pSi 




0 . 0 0 


1 1 584 . 9 


94.5 


0.-25 


9 66 . 1 


1279. 9 


0 . 50 


7335.6 


1501 . 9 


0 . 75 


6 3 4 7 . 


1731 . 7 


1 . 00 


•5655 « 0 


1 7 1 3 . 2 


1 . 25 


5184, 1 


1 703 . 5 


1 . 50 


4832.4 


1 6 26.7 


1 . 75 


4564.7 


1 156 0 . 6 


£ • 0 0 


4361.1 j 


1 1434.7 


XIII. SHALLOW NOTCH FEA LINEAR RESULTS - GAUSS 




a = 4225 psi 
n 




DISTANCE 






FROM NOTCH, in. 


/Psi 


/Psi 


0 . 00 


1 1 530 . 2 


1 1 . '5 


0 . 25 


6 i' i ' • .[ 


i 1 i' i' « r' 


0 . 50 


1 -J J . 


1 504 . 8 


0. 75 


6280 . 9 


1 313. 6 


1 . 00 


5661 . 2 


1 720 . 4 


1 . 25 


5159.3 


1 o29 . 7 


1 . 50 


4 y 4 . 6 


1634. 2 


1 . 75 


4552.8 


1561.7 


2 . 0 0 


4353.6 


1486. 0 
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TABLE XIV, DEEP NOTCH FEA LINEAR RESULTS - NODAL 



a = 4800 psi 
n ^ 

DISTANCE 




TABLE XV. DEEP NOTCH FEA LINEAR RESULTS - GAUSS 

a = 4800 psi 
n ^ 

DISTANCE 
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TABLE XVI. SHALLOW NOTCH FEA NONLINEAR 60,000 lb LOAD 



NO LOAD RESIDUALS 



FROM NOTCH, in. 0 ,psi 


,psi 


»psi 




0 . 0 0 


80024. 7 


1 0 9 4 . 5 


-3341 1 . 4 


439. 1 


0 . 25 


•81338.7 


9523 . 3 


-4250 0 


-2690. 1 


0 • 5 0 


81383.8 


18052 . 8 


9 0 0 8 . 8 


2 1 88 . 3 


0, 75 


b 5 b 4 8 . 8 


1 8993 . 2 


2948. 8 


2133.4 


1 . 00 


58431 . 4 


1 8 799. 8 


1992. 8 


1 798, 7 


1.25 


52870. 1 


1 8 8 47. 


1187,, 3 


1501. 8 


1 . 50 


49111.1 


17477. 1 




1250. 8 


1 . 75 


48080. 7 


18587. 3 


Q ill 


1048.7 


2 , 0 0 


4 3 9 7 3 . 8 


15887,3 


471. 8 


885,0 



TABLE XVII. SHALLOW NOTCH FEA NONLINEAR 65,000 lb LOAD 



DISTANCE 

FROM NOTCH, in. a ,psi 


/Psi 


NO LOAD : 
<^Q >psi 


RESIDUALS 

^Psi 


0 . 0 0 


80390 . 4 


1004.8 


-42092.. i: 


j 224.7 


0.25 


8 2 9 3 9 . 3 


10229. 5 


-9798 . 3 


— 2 I '.3 2 . U 


0 . 5 0 


83842.7 


18054. 8 


8944. 3 


-711.2 


0.75 


74822 . 0 


21475.5 


7507. 3 


3214.2 


1 . 00 


84770. 4 


21171.3 


3 :3 2 0 . ' 


,■'8:38 . 7 


1 . 25 


5789 1 . 2 


20593 . 8 


2202 . 8 


24 1 9. 0 


1 . 50 


5 3 7 8 8 . 7 


19537.0 


1524.5 


2018. 5 


1 . 75 


50321 . 1 


18474.8 


1089, 5 


1710. 1 


2 . 0 0 


47948.2 


17421.4 


1^1 cr j 

‘ .J * a 1 


1448. 0 



37 



TABLE XVIII 



SHALLOW NOTCH FEA NONLINEAR 70,000lb LOAD 



NO LOAD RESIDUALS 



DISTANCE 

FPOM NOTCH, in. ,psi 


0^ rPSi 


fPsi 


fPsi 


0 « 0 0 


01957. '? 


704.5 


~ 5 0 6 7 0 . 6 


6 3 . 9 


0 . ,!1 5 


02822. b 


1 U 6 6 8 . i3 


- 1 6 1 7 6 . !:■ 


-3285. 1 


0 . '5 0 


84007. 0 


15514.4 


-376. 9 


— 2 -i! 6 6 . u' 


0. 75 


0 i‘ 0 1 r . 4 


ii! fci • '9 


15132.6 


2 r' ill. • •-» 


1 . 00 


7 1 2: 1 8 . 0 


237S.J . 2 


5681 . 1 


3983 . 0 


1 . 25 


68782. 9 


23215.2 


■ j 7 6 • j . 1 


b 1 b . b 


1 . 50 


58941 . 5 


21920. 1 


- 6 6 9 n 9 


3030 . 0 


1 . 75 


54897.6 


2 0 670 . 6 


1865,. 4 


259b. 7 


2.00 


52177. 0 


19406. 6 


1 448 . 2 


2207. 0 



TABLE XIX. DEEP NOTCH FEA NONLINEAR 30,000 lb LOAD 



NO LOAD RESIDUALS 



DISTANCE 

FROM NOTCH, in. a ,psi 

9 


a /psi 
r 


fPsi 


a ,psi 
r 


y a y 


8352 1 . 0 


2332.4 


-14138.0 


_ “7 • j cr “i 


0 a 25 


58038 . 6 


20971. 3 


5630. 1 


1402.5 


0 a 5 0 


43309. 1 


18640. 5 


ri IT '7 

•2' '2» -J II r 


i3 24,6 


0^ a 7 5 


34402. 6 


17339. 1 


198. 3 


‘4 6 8 « 9 


1 . 00 


30522. 9 


1 6 0 7 6 . 1 


125.0 


2i 01 01 . 5 


1 . 25 


27543.7 


14574.8 


88 '=1 


215. 2i 


1 . 50 


25757. 0 


1 3328 . 7 


67 . 3 


160. 7 


1 . 75 


24138.8 


12195.0 


5 -I 5 


126.3 


2 . 0 0 


23141 . 6 


11265.7 


4 4 a 9 


100. 9 
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TABLE XX 



DEEP NOTCH FEA NONLINEAR 35,000 lb LOAD 



NO LOAD RESIDUALS 

DISTANCE 

FROM NOTCH, in. ,psi ,psi Cg zpsi /psi 



0.00 


79020 . 4 


1253. 1 




-1909. 


0. 25 


r 'z' cL 1 5 . i' 


24598, 9 


1 fz* z! 4 c* « 7 


1 1 1 :8 . 


0 . 5 0 


50094 . 6 


22675. 1 


1 05 . 7 


1699. 


0. 75 


403 1 b . 0 


21359.5 


624. 1 


1002. 


1 . 00 


35S43. 1 


19156.9 


■3 ci y - 6 


716, 


1.25 


32290. 3 


17299, 1 


257 3 3 


521 . 


1 . 50 


30172.9 


15782. 1 


203. 5 


405. 


1.75 


283 14.8 


14410.8 


1 5 6 - 9 


3 1 9 . 


3 y 0 


27079. 2 


13294.4 


1 2'' 21 1, r' 


261. 



TABLE XXI. 



DEEP NOTCH FEA NONLINEAR 40,000 lb LOAD 



NO LOAD RESIDUALS 



UlbTiUNUi 

FROM NOTCH, in 


CTg ,pSi 


cfj. /Psi 


,PSi 


,psi 


0 3 0 0 


80929 . 5 


835.4 


’ -52900.9 


— 2 1 2 1 . 6 


0 3 25 


89730, 7 


24933 . 8 


20503. 1 


-2132. 3 


0 3 50 


56572. 6 


26815.4 


- 6 1 . 4 


2977 . 0 


0 3 7 5 


4 1 9 84. 6 


24878. 1 


1468. 2 


1643. 1 


1 3 00 


41400. 6 


222 : 0 b . 2! 


6 1 0 . 0 


1 2 1 6 . 6 


1 . 25 


3 7164. 5 


2 0 0 6 3 3 :3 


555.9 


892. 7 


1 . 50 


2' 4 b 9 01 , 9 


13237.2 


438 . 5 


710.5 


1 . 75 


32512.6 


13672- 2 


330 . 7 


564.8 


2 . 00 


31078. 3 


1 5363 . 3 


282. 0 


470. 8 
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mm 



TABLE XXII 



RIGID - PERFECTLY - PLASTIC RESULTS 



a = 73,000 psi 



DISTANCE 



FROM NOTCH, in. 


Og /PSi 


a ,psi 
r 


0 . 00 


73334.4 


J 1 . J 


0. '25 


81730.5 


20363 . 1- 


0 . . J 0 


33252. 6 


2':'' 8 1 3 . ID 


0. 75 


u’ 4 0 r 4 . 1 


37321 . 1 


1 . 00 


34103. 1 


40244. 0 


1 ^ 


3409 1 . 3 


42244 0 


1 . 50 


34051 . 5 


43249.. 3 


1 . 75 


34005 . 0 


'“!• 4 U -.I' ill! II 


2 . 0 0 


33953. 1 


■4 '.:I o 7 ‘li .. *::> 



TABLE XXIII. EXPERIlffiNTAL DATA \ = 0.25 HOLE LINEAR LOADING 

= 10,749 psi 



DISTANCE 



FROM HOLE, in. 


Oq /Psi 


a /Psi 
r 


0 . 000 


35221 . 0 


4 9 . 5 


0.125 


24639 . 0 


- 1650.. 5 


0 . 25 0 


1 9 3 01 9 . 5 


6 6 9 . 5 


0 . -3 7 5 


1 7799. y 


2311 . 5 


0 . 500 


16533.0 


3422.0 
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TABLE XXIV. EXPERIMENTAL DATA SHALLOW NOTCH LINEAR LOADING 

15,000 lb Load 

DISTANCE 

FROM HOLE, in. a ,psi 0 ,psi 

2 T 



0 . 000 


iii! b b 1 1 , 7 


!”i “j ,, V 


0. 125 


22295- 0 




0 . 25 0 


[ q q q ^ j 


• j •”! c." •”( -y 

ii '.j - j c. II r 


01 . 3 7 5 


17277- 5 


2 c! 7 - b 


0 . 5 0 0 


1579S- 7 


2830- 5 


0. 625 


1 ’5 b 7 9 . 1 


1929- 3 



TABLE XXV. EXPERIMENTAL DATA SHALLOW NOTCH 60,000 lb LOAD 



NO LOAD RESIDUALS 



DISTANCE 

FROM NOTCH, in. a ,psi 


,psi 


fPSi 


,psi 


0 . 000 


S 1 1 1 1 . 5 


5551 . 0 


2 3! o 4 2 .. 0 


4.0 


0.125 


79710. 7 


6 7 2 b . b 


~ 1 9 y j 7 , y 


4 1 5 - y 


0 . 250 


73112.0 


3 141.5 


- 1 6409 . 0 


-1044- 0 


0 . 375 


77199.3 


11399.5 


-12620, 01 


-1428.0 


01 . 5 0 0 


r 1 w U r . 2 


1 0509 . 4 


- 1 2 3: 1 9 , ! ;) 


- J 0 0 8 . 0 


01 . 625 


7 1 4 6 . 4 


1 9260 . 0 


- 1 0 9 0 3 . 0 


“ irj 4 5 3 - 0 


0. 750 


63045. 8 


16530. 9 






0 . S75 


59477.4 


1 7 7 4 o . 6 






1 . 000 


57146. 4 


19275.4 
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TABLE XXVI 



EXPERIMENTAL DATA SHALLOW NOTCH 65,000 lb LOAD 



NO LOAD RESIDUALS 



DISTANCE 

FROM NOTCH, in. Q 


»psi 


/psi 


/Psi 


0 . 000 


83073. 9 


7058 . 0 


■ b i '' 4 b . Li 


- 1 '? 1 . 0 


0 . 125 


80957 . b 


6078. 9 


-32543.0 


■•-:j37 . 0 


0 . 250 


79854 . 9 


7468. 1 


-•"23256 . 0 


-1554.0 


0 . 3 7 5 


79482. 8 


10878. 9 


-23914,, 0 


- 233 i , 0 


0 . 5 0 0 


79b738 . 4 


1 6 6 3’ 4 . i j 


-20 1 92 . 0 


-2686 . 0 


U . b 5 


77887. 0 


1 9867 . 6 


- 1 6 3 0 0 . ij 


-3392.0 


0 . 750 


72556.4 j 


19706, 7 






0. S75 


6 4 cJ b p '=b 


17433. 1 






1 . 000 


b' 1 4 c 5 y ^ 


: 19034.1 







TABLE XXVII. EXPERIMENTAL DATA SHALLOW NOTCH 70,000 lb LOAD 

NO LOAD RESIDUALS 



DISTANCE 



FROM NOTCH, in. 0 . ,psi 

o 


»psi 


Og ,psi 


,psi 


0 . 000 


84743.9 


8 569. 0 


- 5 0 7 9 1 . 0 


43.0 


0. 125 


93492.9 


9073. 9 


-46086 . 0 


-736. 0 


0 . 25 0 


84258 . 3 


1 

137 18. 2 


1 

-40774.0 


— c* r :!•!: . 01 


0 . 375 


92846.8 


1 3:8 1 5 . 0 


- :2 4 6 0 0 , 0 


- 1 0 1 15 , 0 


01 . 5 0 01 


92524. 6 


1 6 :5 1 4 , 9 


- 3 0 0 1 7 . 0 


- 1 2 : 1 16 . 0 


0. 625 


9 0446. 4 


17441 . 6 


-24798 . 0 


- 1 551. 0 


0. 750 


90455.5 


23190.4 






0 . 6* i' 5 


76569 . 9 


229 2:'' . 7 






1 . 0010 


72264 . 1 


23604 . 0 
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TABLE XXVIII 



EXPERIMENTAL DATA DEEP NOTCH LINEAR LOADING 



15,000 lb Load 



DISTANCE 



FROM NOTCH, in. 


rPSi 


a fpsi 
r 


0 . 080 


45 '? 07. 0 


- 45. 0 


0 . 185 


35372 . 0 


7693 . 0 


0 . 250 


3 0 410. 0 


lit 7 1 „ 0 


0. 375 


2 ' 6 1 7 7 . 0 


1 2395. 0 


0 . 500 


20933 . 0 


9879. 0 


0. 625 


13373. 0 


9060 . 0 


0. 750 


1 689 1 . 0 


3934. 0 



TABLE XXIX. EXPERIMENTAL DATA DEEP NOTCH 30,000 lb LOAD 

Elastic - Plastic 

DISTANCE 



FROM NOTCH, in. 


Vg ,psi 


a ,psi 
r 


0.000 


76156. 6 


• J '2' i -I- 


0.125 


00996.2 


6 0 57.7 


0 . 250 


50047. 3 


[ 0456. 3 


0. 375 


4 0 9 9 2 . 0 


1 4435 . 3 


0 . 500 


2 ' 9 7 4 cl . 6 


1643 3 . 6 


0 . 625 


35030 . 8 


1 5399 . 21 


0 . 750 


cl cl 2 0 6 . 0 


16720. 0 
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APPENDIX A 



PSAPl JCL 



STANDARD J3B CARO 

// EXEC FRTXCL3P 

//FORT.SYSPRINT DO DUMMY 

//FORT.SYSIN DO UN IT=3330 ,VOL=SER=0 I SK02 t 

// DSN=S2 939,PSAP(PSAP) . 0 IS P=SHR» LA 8EL= ( » » » IN) 



// 


DO 


UN 


IT= 


3330 


♦ VOL 


=SER = 


0 


// 


DISP= 


SHR 


♦ LAB 


EL=( 


♦ ♦♦IN 


) 


// 


OD 


UN 


IT= 


3330 


♦ VOL 


=SER = 


0 


// 


OISP= 


SHR 


♦ LAB 


EL = ( 


♦ ♦♦IN 


) 


// 


DO 


UN 


IT= 


3330 


♦ VOL 


= SER = 


0 


// 


01 SP= 


SHR 


♦ LA3 


EL = ( 


♦ ♦♦IN 


) 


// 


OD 


UN 


IT= 


3330 


♦ VOL 


=SER = 


0 


// 


OISP= 


SHR 


♦ LAB 


EL=( 


♦ ♦♦IN 


) 


// 


DO 


UN 


IT= 


3330 


♦ VOL 


=SER = 


D 


// 


OISP= 


SHR 


♦ LAB 


EL = ( 


♦ ♦♦IN 


) 


// 


DO 


UN 


IT= 


3330 


♦ VOL 


=SER = 


D 


// 


01 5 


;p= 


SHR 


♦ LAB 


EL=( 


♦ ♦♦IN 


) 


// 


DO 


UN 


IT= 


3330 


♦ VOL 


=SER = 


0 


// 


OISP= 


SHR 


♦ LAB 


EL = ( 


♦ ♦♦IN 


) 


// 


DO 


4c 












C 


MA 


IN 


PROGRAM 


4c4c4c 








DI 


MENSION 


zzz 


(3203 


) 






CALL 


PSAP 


1( zz 


Z^320 


3 



STOP 

END 



SK02t0SN=S2939. PSAP (PLOT) 
SK02tDSN=S2939.PSAP( IN IT) 
SK02fDSN=S2939.PSAP(ELER) 
SK02,DSN=S2939. PSAP( SAPF) 
SK02tOSN=S2939.PSAP( ADNA) 
SK02,0SN=S2939.PSAP( AUXL) 
SK02,0SN=S293 9.PSAP( AOPT) 

OISPD(5,3,80O) 

OISPO,3O0 ) 



DELIMITER CARD (/*) 

//GD.FTlOFOOl OD UN IT = SYSDA, 0 ISP=( » PASS ) t 
// SPACE=(CYL,( 2»2) ) ,DSN=SCTEMPl 
//GO. SYS IN DO * 

:4c « 4 4c 3jc 4c 4c # # 4c ^ « :(c:(c:{c « 4c 4c 4c 4c 3)c 4c 4c # ^ ^ 30c :4c :4c 4: :(c :4c 4c 4e :(c 

4c4c 4c 4c 4c 4c 4c 4c 4c INSERT PSA^l DATA HERE 4c 4c 4c 4c 4c 4c 4c 4c 4c 4t 4c 

4c4c4c4c4c4c4c4c4c4c4c4:4c4cS4c4c4c4c 4c4c4c4c4: 4c 4c4c4c4c 4c4c4c4c 4c4c 4c4c4c4c 4c4c4c4c4c4c 
4c4c4c4c4c4c4c4c4c4c4c4c4c DELIMITER CARO (/*) 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4= 4c 4c 4c 



f 
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APPENDIX B 



LOCAL LEAST SQUARES SMOOTHING 




• Corner Nodes 
o Midside Nodes 
+ 2x2 Gauss Points 



Two-Dimensional Isoparametric Element from ADINA [Ref. 4] 



The local smoothing expression from Hinton and Campbell 
[Ref. 19] in ADINA coordinates becomes 
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With 02/ 0^ 3.nd 5^ representing the smoothed corner node 

stresses and 0 ^, -[•-[•/ the unsmoothed 

stresses at the Gauss integration points, this expression 
can be written in an equivalent form. 
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The rnidside node stress values may be obtained by averaging 
the values at the associated corner nodes, since the distri- 
bution of the smoothed stresses is linear along the sides of 
the element. Smoothed stress values obtained by this least 
squares method should subsequently be averaged to obtain 
unique values at nodal points shared by adjacent elements. 
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APPENDIX C 



ADINA JCL 



STAMOARO JD3 CARO 

// EXEC FCRTXCLG,RE3I3N = 200-:)K 
//FORT.SYSPRI^J^ DO OJ'I'iy 
//FORT.SYSINI DO * 

IMPLICIT REAL*8( A-H,0-Z ) 

REAL A 

COMMON A( 120010) 

MAX=1200D0 
CALL EXE: (MAX ) 

STQD 

END 

^ ^ ^ ^ DELIMITER CARO ( /* ) ^**^**^^^f*-^**^^* 

//LKED.JSOO DD DI SP = SHR , D SiN= M SS . S 2 9 39 . AOI NA 
//LKED.SYSIN DD * 

INCLUDE 'JSOD(LOADM) 

ENTRY MAIN 

DELIMITER CARD { /^) 

//GO. FT07F001 DO UN IT = S Y SOA, D ISP= ( N EtJ , 0 ELET>= ) , 

// OCB=(RECFM=VBS »3L<SI ZE=A0D0) , S P ACE = ( CYL r ( 5 , 1) ) 
//GO. FTOlFOOl DD UN I T = S Y SOA , 0 ISP = ( N EW , 0 ELETE ) , 

// 0CB=(RECFM=V3SfBL:<5IZE=A3 0D) ,SPACE=(CYLt (5,1)) 
//GO. FTD2F001 DO UN I T = S '' SDA , D ISP= ( N EW , 0 ELETE ) , 

// OCB=( RECFM=VBS,BL<5IZE=^D00) , S? ACE = ( C YL , ( 5,1)) 
//GO. FTD3F001 DO UN IT = S YS OA , 0 ISP= ( N EW , D EL ETE ) , 

// 0C8=(RECFM=V3S,3L<SIZE=1D00) ,S?ACE=(CYL,( 5, 1) ) 
//G0.FT04F001 DD UN IT = S Y SDA , D ISP = ( N EW , DE LETE ) , 

// 0C8= ( RECFM= V8S,BL.<SI ZE = 1DD0) ,SPACS= (CYL, ( 5 , 1 ) ) 

//GO. FT08F001 DO UN IT = S / SOA, D ISP= ( N EW , 0 ELE TE ) , 

// 0C3=(RECFM=V3S,3L<5IZE=A3D0) , SPACS= ( CYL, ( 5 , 1 ) ) 
//GO. FT39F001 DD UN I T = S Y SOA , 0 IS? = ( NEW , DE LETE ) , 

// 0CB=(RECFM=VBS,BL<SIZE=^0 00),SPACE = (CYL,( 5,1)) 
//GO. FTlOFOOl DD UN IT = S Y SOA , 0 ISP= ( N EW , 0 ELETE ) , 

// 0C3=(RECFM=VBS,3LKSIZE='+0 00) , SPACE= ( CYL, ( 5,1)) 
//GO.FTllFOOl DO UN IT = S Y SOA, D ISP= ( N EW , 0 ELETE ) , 

// 0C3=(R5CFM=V3S ,3LKSI ZE = 13 00) , S P ACE = ( C YL , ( 5 , 1 ) ) 
//GO. FT12F001 DD UN IT = S Y SOA , D ISP= ( N EW , DE LETE ) . 

// OC3=(RECFM=VBS,3L<SIZE=13DO),S?ACE=(CYL,(5, 1) ) 
//GO. FT13F001 DO UN I T = S Y SOA , 0 ISP= ( N EW , DELETE ) , 

// 0C8= ( RECFM=V8S, BL<SI ZE=1000) , S?AC£ = (CYL, (5,1)) 

//GO. FT56F001 DO UN IT = 3 3 30 , VC L=SER= 01 SKO 1 , 

// 0SN=F0 399. TEMP,DISP=SHR,LA8EL=( , , , I N) 

//GO. FT57P001 00 UN IT = SY SOA, D ISP= ( NEW , DELETE ) , 

// 0C3=( RECFM= VBS,3L<SIZE=13D0) , SPACE= ( CYL, ( 5 , 1 ) ) 

//GO.SYSIN 00 

♦ :(C *3je # * * 

INSERT AOINA DATA HERE 

BLAN< CARD 

BLANK CARO (TWO BLANK CAROS STOP EXEC) 

DELIMITER CARD (/*) 
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APPENDIX D 
PSAPl LISTING 
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